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We develop the multi-scale ab-initio scheme for correlated electrons (MACE) for transition- 
metal-oxide heterostructures, and determine the parameters of the low-energy effective model. 
By separating Ti t2g bands near the Fermi level from the global Kohn-Sham (KS) bands 
of LaAlOs/SrTiOs which are highly entangled with each other, we are able to calculate the 
parameters of the low-energy effective model of the interface with the help of constrained 
random phase approximation (cRPA). The on-site energies of the Ti t2g orbitals in the 1st- 
layer is about 650 meV lower than those in the 2nd-layer. In the Ist-layer, the transfer integral 
of the Ti t2g orbital is nearly the same as that of the bulk SrTiOa, while the effective screened 
Coulomb interaction becomes about 10 percent larger than that of the bulk SrTiOs. The 
differences of the parameters from the bulk SrTiOs reduce rapidly with increasing distance 
from the interface. Our present versatile method makes it possible to derive effective ab-initio 
low-energy models and allows studying interfaces of strongly correlated electron systems from 
first principles. 

KEYWORDS: first-principles calculation, effective Hamiltonian, downfolding, constrained RPA method, 
correlated-electron systems, heterostructure, interface, two-dimensional electron systems 
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1. Introduction 

In recent years, interfaces of strongly correlated 
electrons have been under intense investigations. Es- 
pecially, transition-metal-oxide heterostructure SrTiOa 
(ST0)/LaA103 (LAO) has received a lot of attention, 
because of its remarkable transport properties.^ The in- 
terface of the SrTiOa/LaAlOs shows metalhc conductiv- 
ity, although the bulk materials of each transition-metal- 
oxide, SrTiOa and LaAlOs, are band insulators. The 
Ti02-terminated (n-type) interfaces show the metallic 
conductivity, when the thickness of deposited LAO layers 
is thicker than four unit cells, ^ but the SrO-tcrminatcd 
(p-typc) interfaces are insulating for any LAO thickness.'^ 
The origin of this conductivity at the n-type interface 
remains a debated issue. An intrinsic electronic effect, 
namely the polar discontinuity, and an extrinsic atomic 
effect, atomic vacancies, were proposed to play key roles 
at high carrier concentrations. When the polar LAO lay- 
ers are deposited on the Ti02-terminated substrate STO, 
an electric potential along the [001] direction diverges as 
the LAO thickness increases. The instability of this elec- 
tric potential is suppressed, if 1/2 doped electron or a 
corresponding atomic charge vacancy per unit cell ex- 
ist at tlie interface.'^'* Besides such transport properties, 
the existence of superconductivity and magnetic order 
are also reported at the SrTiOs/LaAlOs interfaces.^' ^ 
These unique transport properties are expected to offer 
a useful functionality in the possible applications. 

The Local Density Approximation (LDA) calculations 
for transition-metal-oxide heterostructures have already 
been performed by several groups. Park et al. inves- 
tigated LaAlOa/SrTiOs for three kinds of supcrlattice 
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structures with n-type. p-type, and both types of in- 
terfaces.^ By taking into account the relaxation of lat- 
tice, Ishibashi and Tcrakura investigated the influence of 
LaAlOa thickness on the carrier density at the interface.* 
The LDA, however, often fails to capture correlation ef- 
fects in trairsition-metal-oxides. To study novel electronic 
phases in transition-metal-oxide heterostructures, in ad- 
dition to the carrier doping from the electronic recon- 
struction, we should treat the correlation effects beyond 
the LDA, because the correlation effect is, on general 
grounds, expected to be enhanced at interfaces because 
of the effective reduction of the spatial dimensionality. 

The LDA-f U is useful for the strongly correlated ma- 
terials, especially for insulators. Pentcheva and Pickett 
tried to explain the insulating behavior of the p-type 
LAO/STO interface by introducing the on-site U, since 
the LDA gives a metallic state. ^ However, in the LDA-I-U 
method, there is no established ab-initio way to estimate 
U, although the value of U strongly affects the calculated 
results of the low-energy physics. Furthermore, owing to 
the experimental difficulties, U has to be cited from the 
experimental values of the bulk, not of the interface, es- 
pecially in the case where thick LAO layers are on a STO 
substrate. Furthermore, the single-particle theory often 
collapses in strongly correlated materials, even with the 
suitable U. 

To overcome such serious problems, a combined 
method of the multi-scale ab-initio scheme for corre- 
lated electrons (MACE) is very useful.^" With the con- 
strained random phase approximation (cRPA), we can 
calculate not only the transfer integral but also the 
screened Coulomb interaction without relying on any ex- 
perimental parameters. Correlation effects can be treated 
accurately by solving the obtained ab-initio effective 
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model using low-cncrgy solvers such as variational Monte 
Carlo, dynamical mean field theory (DMFT)^^'^^ 
and path-integral renormalization group (PIRG).^'^' 
The MACE has already been applied to a wide vari- 
ety of materials; semiconductor/^ transition-metal,^''"^^ 
transition-metal-oxide, ^^"^'^ molecular organic conduc- 
tors^'* and Fe-based layered superconductor.^^' So far, 
the MACE has been applied to the bulk material. Since 
interfaces have opened a new avenue of research, it is 
highly desired to develop a methodology which is able to 
treat the interface from the same footing of the MACE to 
understand the electron correlation effect from the first 
principles. 

In this paper, we apply the ab-initio downfolding 
method to transition-metal-oxide heterostructures, and 
determine an ab-initio low-energy effective Hamiltonian 
of the SrTiOa/LaAlOs. We derive the effective model 
from the conduction bands near the Fermi level origi- 
nated mainly from the Ti t2g orbital. The on-site energies 
of the Ti t2g orbitals in the Ist-layer is about 650 meV 
lower than those in the 2nd-layer. In the Ist-laycr, the 
transfer integral is nearly the same as that of the bulk 
SrTiOs; while the effective screened Coulomb interaction 
becomes about 10 percent larger than that of the bulk 
SrTiOa. The parameters of the 2nd-layer are similar to 
those of the bulk SrTiOa. 

In Sec. 2 we describe our method. Section 3 describes 
the band structure and the derived effective model of the 
SrTiOs/LaAlOs. We also present an effective model of 
SrTiOa for comparison. Section 4 is devoted to summary. 

2. Methods 

To construct the effective low-energy model of the 
heterostructure with well defined parameters, we treat 
the screening by the high-energy parts working on the 
low-energy parts starting from the band structure of 
LAO/STO, and renormalize the high-energy parts into 
the low-energy parts. This downfolding procedure is in- 
troduced by Aryasetiawan et al?"^ and Solovyev et al}'' 
The low-energy effective model obtained from the down- 
folding procedure offers a starting point for studies on 
low-energy physics. 

We consider an extended Hubbard Hamiltonian de- 
scribing low-energy electronic properties of the interface 
of LAO/STO, 

^ = X] X! X] ^ Rj)"'il"-'jm 

cr ij nm 

ap ij nm ^ 

+ Jrnn{Ri ~ -Rj ) K°2a^maf„ajm + «rj«m«im'^jm) | > 

(1) 

where afj^ (af„) is a creation (annihilation) operator of 
an electron with spin a in the n-th orbital, which is de- 
fined by a maximally localized Wannier function centered 
at the ili-th unit cell.^®'^^ Especially, in this paper, we 
construct a model for all Ti t2g orbitals, not only at the 



interface but also in the bulk STO region, in the super- 
cell of LAO/STO. With thick LAO and STO layers, this 
3-dimcnsional model corresponds to the semi-infinite 3- 
dimcnsional model for (vacuum)-(LaA103)Ar-(SrTi03)oo 
where the LAO layers are deposited on the substrate 
bulk STO, because the inter-supercell screening effect is 
rather weaker than the intra-one. Therefore, the param- 
eters in this calculation may be used not only for super- 
lattices but also for a semi-infinite interface LAO/STO. 
The parameters of the Ti t2g orbitals away from the in- 
terface recover the bulk STO nature in both the super- 
cell and semi-infinite interface, because the polarizations 
away from the bulk STO region contribute little to the 
screening to such t2g orbitals. In fact, in the present cal- 
culation, we show that the parameters in the 2nd STO 
layer from the interface nearly converge to those of the 
bulk STO. With thick LAO layers, the parameters of 
this 3-dimensional model are also similar to those of a 2- 
dimensional model for LAO/STO superlattice, where all 
of the inter-supercell screening effects are renormalized. 

We derive the effective model in the following way: 
First, we calculate the whole band structure of the 
transition-metal-oxidc heterostructure in the framework 
of density functional theory (DFT). Wc then choose 
the target low-energy band around the Fermi level and 
construct the maximally-localized Wannier functions^^ 
in the low-energy Hilbcrt space. The transfer integral, 
which defines the one-body part of the low-energy model, 
is obtained as the matrix element of the Kohn-Sham (KS) 
Hamiltonian Hks, 

tmn{R) = ((^Oml'HKsl^fin), (2) 

where 0ij„ is the n-th Wannier function centered at the 
ceU R. 

Next, we renormalize the effect from the high-energy 
space (r space) into the target low-energy space (d 
space), and evaluate effective electron interaction by us- 
ing the cRPA. We divide the polarization 

occ unocc 

P(r,r» = ^ ^ ^,(r)^*(r')^/'*(r)^,(r') 

i 



-U) — tj + ei + i5 uj + ej — ei — i5\ 

into the polarization Pd that includes only the d-d tran- 
sitions and the rest of the polarization ■ The screened 
Coulomb interaction W is given by 

W = [l-vP]-'^v 

(4) 

= [\-WrPdVWr, 

where we define the partially screened Coulomb interac- 
tion Wr that does not include the screening arising from 
the polarization from the d-d transitions as 

Wr^[l-vPr]~^V. (5) 

The screened Coulomb matrix is defined by 
Wr{R\n, R-^n' , Rzm, R^m'; w) 
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= j dVdV'0l,^„(r)</>K,„'(r) 

X W^,(r,r';c.)0^3„(r')<^H.™'(r'). (6) 

The effective Coulomb interaction U and the exchange 
interaction J are given by 

UnniiR) = WriOn,On, Rm, Rm;0), (7) 

Jnm{R) = Wr(On, Om, Rm, Rn; 0). (8) 

In LaAlOa/SrTiOs, the Ti states are entangled 
with the La 4/ states. We disentangle the Ti bands 
using the recently developed disentangling technique. 

In the LDA, the La 4/ level in LAO/STO is located 
near the Fermi level. However, in the real material, the 
La 4/ bands are expected to be at higher energy, and, as 
a consequence, do not strongly screen the interactions be- 
tween the electrons in the t2g bands. To calculate the pa- 
rameters of the low-energy model with a higher accuracy, 
a better way is to take into account the correlation effect 
of the La 4/ bands by the GW approximation (GWA) as 
a preconditioning before determining the effective model 
for Ti M t2g orbitals. The detail of the GWA for the 
disentangled 4/ bands is explained in Appendix. 

Computational conditions are as follows. We calcu- 
late the band structures of the bulk STO and LAO, 
and the LAO/STO heterostructure based on the DFT- 
LDA.^^'^^ The calculations are carried out with the pro- 
gram based on the full-potential linear muffin-tin orbitals 
(FP-LMTO) method. 3'"' Because the LAO are deposited 
to fit the substrate STO, the lattice constants of the de- 
posited LAO layers in the direction horizontal to the in- 
terface are same as those of the substrate STO, while that 
in the direction perpendicular to the interface slightly 
changes depending on the carrier density at the inter- 
face. Therefore, in this study, the structure is fixed as 
cubic, and the lattice constants for the bulk and the het- 
erostructure perpendicular to the [001] stacking direction 
are fixed to 3.905A, which corresponds to the experimen- 
tal lattice constant of the bulk SrTiOa. In the LDA cal- 
culations, 8x8x8 fc-point sampling is employed for the 
bulk, and 8x8x2 fc-point sampling is employed for the 
LAO/STO to represent electronic structures of the sys- 
tem. The muffin-tin (MT) radii are: R^^ 2.50 bohr. 
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and i?o"^ = 1.6 bohr. The angular momentum cut off is 
taken at / = 4 for all the sites. In the cRPA and the GW 
calculations, 3x3x3 fc-point sampling is employed for 
the bulk, and 3x3x1 fc-point sampling is employed for 
the LAO/STO. 

3. Results 

3. 1 band structure and density of state 

First, we show the band structure and the density of 
states of the insulators, bulk SrTiOs and LaAlOa. The 
upper panels of Figs. 1 and 2 show the band structure 
and the density of states , respectively. At room temper- 
ature, SrTiOs has the cubic perovskite structure, and 
becomes tetragonal below 105 K. SrTiOs has a high di- 
electric constant at low temperatures because of the na- 
ture of quantum paraelectricity.'^^ In this calculation, the 
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Fig. 2. (color online) Density of states of bulk STO and LAO, 
and LAO/STO obtained by LDA. Energy is measured from the 
Fermi level. 



structure is fixed as cubic, and the lattice parameters 
are fixed at 3.905 A of the bulk SrTiOg. The lower three 
conduction bands are derived from the orbital of Ti 
sites, where the octahedral crystal field of 0^~ partially 
breaks the 5-fold symmetry of the 3d orbitals into the 
lower orbitals of the t2g and the higher orbital eg. The 
calculated band gap is 2.1 eV and the band width of 
the t2g band is 2.8 eV. In the experiment, the band gap 
is 3.3 eV.^^ In the strongly correlated materials, such 
underestimation of the band gap causes to overestimate 
the screening effect from the high-energy bands to the 
low-energy bands in the cRPA. In the LAO/STO, the 
underestimation of the energy levels of La 4/ bands is a 
major problem which are entangled with Ti t2g bands in 
the LDA as we will show. LaAlOa has the rhombohedral 
perovskite structure at room temperature, and becomes 
cubic above 821 K. To compare with the LAO layer in the 
heterostructure, the structure and the lattice parameters 
are fixed at those of the cubic SrTiOs in this calculation. 
The lower narrow conduction bands of LAO at ^4 eV are 
derived from the 4/ orbital of La sites. The La 5d bands 
are hybridized with the La 4/ bands at the F-point. In 
the experiment, the band gap between the La 4/ and the 
O 2p is 5.6 eV.'^^ In this calculation, however, the band 
gap is 3.4 eV because the LDA underestimates the value 
of the band gap. This band-gap problem in the LDA is 




Fig. 1. Electronic band structures obtained by the LDA. The zero energy corresponds to the Fermi level. Upper left panel: Electronic 
band of STO. Upper right panel: Electronic band of LAO. Lower panels: Electronic bands of LA01.5ST03.5. The right one is the 
enlarged view around the Fermi level. Wave functions of the states indicated by arrows are displayed in Fig. 4. 
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Fig. 3. (color online) Schematic picture of the model of the 
transition-metal-oxide heterostructure LA01.5ST03.5. The sys- 
tem is modeled by a periodically repeated supercell containing 
two interfaces. The both interfaces are LaO-Ti02. Owing to neu- 
trality of the electrons and atoms, these systems are formally 
charged by — e electron per supercell. 



improved dramatically with the GWA. We will show a 
result in the GWA later. 

Next, we show the LDA results of the transition-metal- 
oxide heterostructure LA01.5ST03.5. We refer to the 
Ti02 layer at the interface as the Ist-layer and the Ti02 
layer in the bulk region of the SrTiOa as the 2nd-layer. In 



this paper, we refer to the heterostructure -(LaAlOa)!- 
LaO/Ti02-(SrTi03)3- as LA01.5ST03.5 (see Fig. 3). 
This heterostructure has two crystallographically equiv- 
alent n-type interface, and has about 1/2 carrier electron 
at each n-type interface, because LA01.5+^ST03.5° has 
a positive charge in the ionic limit. In terms of the po- 
lar discontinuity, these models of the heterostructures 
are the cases where the instabilities of the potential di- 
vergence are completely suppressed.'^' ^ The lower pan- 
els of Fig. 1 show the band structure of the n-type 
LA01.5ST03.5 heterostructure. The energy bands are 
rather degenerate because there are two crystallograph- 
ically equivalent n-type interfaces. The LA01.5ST03.5 
has the Fermi surface around the F point. The density of 
states of the LAO/STO comes from nearly superimposed 
states of the bulk LAO and STO (see Fig. 2). 

Figures 4 shows the isosurface contours of selected 
Bloch functions of the conduction bands at the F-point. 
The panel (a) of Fig. 4 is the isosurface contour of the 
lowest conduction band (see Fig. 1 (d)). This band is 
originated mainly from the Ti dxy in the Ist-layer. The 
panel (b) is the isosurface contour of the Ti dyz band, 
which spreads in the direction horizontal to the inter- 
face. The panel (c) and (d) are the isosurface contours 
of the Ti dxy band in the 2nd-layer and the La 4/ band, 
respectively. The Ti <2g and La 4/ orbitals are spatially 
close and, in the LDA level, energetically close, so that 



J. Phys. Soc. Jpn. 

{a)-0.7250 eV 



Full Paper 



Author Name 



{b)-0.1837 eV 




0044 eV 




Fig. 4. (color online) Isosurfacc contours of the Bloch functions of 
LA01.5ST03.5 conduction bands at P-point [000]. The structure 
corresponds to Fig. 3, for example, the panel A represents the 
Ti dxy at the interface. Value of the isosurface is set to ±0.05 
bohr"^/^. Each energy level is measured from the Fermi level 
and is depicted as arrows in Fig. 1. 



the Ti t2g and La 4/ orbitals are hybridized at the in- 
terface (see Fig. 4). The Bloch functions transfer to the 
bulk region with increasing those energy levels. There is 
a positive crystal field from the polar perovskite LaAlOs 
to the non-polar substrate SrTiOa as compared with the 
case with only Sr^"*" for the bulk SrTiOs. This crystal 
field vanishes in the bulk region of SrTiOa, because the 
negative field from the doped electrons compensates this 
positive field. From these reasons, the energy levels of 
the orbitals in the Ist-layer arc lower than that in the 
2nd-layer, and the gap between the valence and conduc- 
tion bands in LA01.5ST03.5 is smaller than that of the 
bulk STO. Such effect of compensation and confinement 
becomes strong as the polar perovskite LaAlOa becomes 
thicker. This tendency is also seen in Fig. 5 which shows 
the partial density of states at the Ti sites of d orbitals 
and the O sites of p orbitals in the Ti02 layer of the 
bulk STO and LAO/STO obtained by the LDA. Here- 
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Fig. 5. (color online) Partial densities of states of Ti d and O p at 
Ti02 layers of bulk STO and LAO/STO obtained by the LDA. 
Energy is measured from the Fermi level. 



after, for instance, a Wannier orbital of the xy orbital in 
the Ist-layer is denoted by Ixy. Pronounced peak shifts 
of DOS to the lower energy are found for the sites in the 
Ist-layer. 

Considering the experimental band gap of the bulk 
LaAlOa, we note that the location of the energy of the 
La 4/ bands in the LaAlOs/SrTiOa is too low in the 
LDA. The La 4/ bands screen and hybridize with the Ti 
t2g bands weaker in the real material. To calculate the 
parameter of the low-energy model with a high accuracy, 
a better way is to take into account the correlation effect 
of La separately by the GWA. To calculate the self ener- 
gies of the La 4/ bands, we first construct the maximally 
localized Wannier function of La 4/ from a linear com- 
bination of the target low-energy KS-bands. We choose 
the energy window from —1 eV to 3 cV for the Wannier 
functions. We find, however, that the screened Coulomb 
interaction of the low-energy bands derived from the Ti 
t2g is not sensitive to the choice of the energy window if 
the window exceeds a certain width but is not too wide. 
Next, we calculate the self-energy corrections of the La 
4/ bands by the 1-shot GW scheme. Figure 6 shows the 
band structure after the 1-shot GW corrections for the 
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Fig. 6. (color online) Electric band structures with 1-shot GW 
self-energy for 4/ bands. Energy is measured from the Fermi 
level. 




Fig. 7. (color online) Disentangled d-bands having strong t2g 
character ((blue) dashed line) and diagonalized r-bands (solid 
line). 



4/ bands. The energy level of 4/ bands is raised by about 
1.5 eV with GW. This value is consistent with the ex- 
perimental result of the LAO.'^^ In the following sections, 
we show this self-energy effect on the effective model pa- 
rameters. 

3.2 Wannier function and transfer integral 

We calculate the transfer integral t of the Ti t2g or- 
bital to determine the 1-body part of low-energy effec- 
tive Hamiltonian for the heterostructure. First, we con- 
struct 6x2 maximally localized Wannier functions hav- 
ing strong Ti t2g characters, from the linear combination 
of the target low-energy KS-bands. We choose —1-3 eV 
as the energy window for the Wannier functions. Figure 
7 shows the disentangled t2g bands and the rest bands, 
and Fig. 8 shows the isosurfaces of Wannier functions of 
the dxy and dyz orbitals. 

We show in Table I the transfer integrals t of the bulk 
STO calculated in the LDA. In the tables and this sub- 
section, m and n specify symmetries of the Ti t2g or- 
bitals. The values of on-site energies are listed in the 
column for {R^, Ry, Rz) = (0,0,0). The values of the 
on-site energies are the same for all the t2g orbitals be- 
cause of the cubic crystal symmetry. The major values 
of the nearest hopping between the same symmetries 




Fig. 8. (color online) Isosurface of the maximally localized Wan- 



nier functions ±0.05 a.u. 
Ist-laycr. 



for the dxv and dyz orbitals in the 



{Rx, Ry, Rz) = (1,0,0) are 299 meV. These values are 
consistent with the values in the literature in nearest- 
neighbor tight-binding models of the transition-metal- 
oxide interfaces {t ^ 0.3 eV).'^^ The nearest hopping for 
the perpendicular directions of the orbitals symmetry is 
—39 meV, which is 13 percent of these of the main direc- 
tions. The next nearest hopping {Rx, Ry, Rz) = (1,1,0) 
is 35 percent of the nearest neighbor hopping. 

Next, we show in Table II the transfer integrals t of 
the LA01.5ST03.5 calculated in the LDA and the GW, 
respectively. The values of on-site energy in the Ist-layer 
iim.im(0, 0, 0) are about 650 meV lower than those in 
the 2nd-laycr t2m.2m(0, 0, 0), mainly because the dipole 
moment of LaO^^-AlOj^ stabilizes the energy of the t2g 
orbitals in the Ist-layer. The hybridization between the 
La 4/ and Ti t2g orbitals also slightly stabilizes the en- 
ergy of the t2g orbitals. These layer dependent potential 
localizes carriers at the interface. Actually, in the exper- 
iment, the transition into the 2D superconducting state, 
namely the Berezinskii-Kosterlitz-Thouless transition is 
seen at the interface of the LAO/STO.^ The value of 
the on-site energy of Ixy, ti2;y.ia:i,(0, 0, 0), is lower than 
that of the other lt2g due to the crystal field and the 
hybridization with the LAO layer. In the GWA, this dif- 
ference of the on-site energy at the Ist-layer is 16 meV 
larger than that in the LDA, because the hybridizations 
of the La 4/ orbitals with the Ti t2g orbitals in the Ist- 
layer, especially with lyz and Izx, become weaker. In 
the 2nd-layer, the dipole moment of LaO^^-AlOj" have 
little effect, and the values of on-site energy partially 
recover the bulk STO nature. The hopping parameters 
are similar to those of the bulk STO. The major val- 
ues of the nearest hopping between the same symmetries 
iim.im(l, 0, 0) are about 0.3 eV. The nearest hopping of 
the Ixy is smaller than that of the bulk STO. The near- 
est hopping of the lyz and Izx are nearly the same as 
those of the bulk STO. The main difference of the hop- 
ping parameters from the bulk STO is seen in the nearest 
hopping for the perpendicular directions of the orbitals 
symmetry tiyz,iyz(l, 0, 0). Such hoppings are —59 meV, 
which are about twice as large as those of the bulk STO. 
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In the LAO/STO, therefore, the carriers tend to spread 
in the direction horizontal to the interface compared to 
the bulk STO. As with the on-site energy, the hopping 
parameters in the 2nd-layer t2m,2m(l, 0, 0) are nearly the 
same as those of the bulk STO. In the GWA, the hy- 
bridization between the La 4/ and t2g orbitals in the 
Ist-layer, especially lyz and Izx, becomes weaker. As a 
result, the splitting of on-site energies of orbitals be- 
comes larger, and the hopping parameter between Izx- 
2zx becomes smaller. 

3.3 screened Coulomb interaction 

Next, we calculate the on-site screened Coulomb in- 
teraction J7„„ (0,0,0) and the on-site screened exchange 
interaction Jmn(0, 0, 0) of the Ti t2g orbital to determine 
the 2-body part of the low-energy effective Hamiltonian 
for the heterostructure. The screened Coulomb interac- 
tion is computed with the matrix elements in the max- 
imally localized Wannier basis in the framework of the 
cRPA (see eq. (6)). 

In the top of Table III, we show the effective on-site 
Coulomb interaction U and J for the bulk STO. As with 
the transfer integrals, m and n denote symmetries of the 
t2g orbitals of Ti in the tables and this subsection. The ef- 
fective on-site Coulomb interaction U between the same 
orbitals is 3.76 eV, while the bare on-site Coulomb in- 
teraction is 14.27 eV. The on-site Coulomb interaction is 
reduced about to 1/4 by the polarizations without the d- 
d contributions. The effective on-site screened exchange 
interaction J is 0.46 eV. 

We show the effective on-site Coulomb interaction U 
and J of LAO/STO calculated with cRPA in the mid- 
dle and bottom of the Table III. Here, U between the 
same orbitals are 3.4-3.6 eV from the LDA and 3.7-4.0 
eV from the GWA, while that of the bulk STO is 3.76 
eV. In both the LDA and GWA, Uixy,ixy is the largest 
among on-site screened Coulomb interactions in the Ist- 
layer, and Uiyz,iyz and Uizx,izx are 4-5 percent smaller 
than Uixy,ixy In the GWA, U and J at the Ist-layer 
become larger than those of the bulk STO, while these 
are smaller in the LDA. This is because the La bands 
are raised away from the Fermi level by the self-energy 
correction, the polarizations between La 4/ and Ti t2g or- 
bitals become weaker compared to the LDA bands, and. 
as a result, the screening effect becomes weaker. The on- 
site screened Coulomb interaction Umn and the on-site 
screened exchange interaction J,„„ satisfy the equation 
Umn = Umm — '^Jmn in thc bulk STO. On the other hand, 
at the interface of the LAO/STO, these parameters do 
not satisfy this equation because of the inversion sym- 
metry breaking. In the 2nd-layer, U and J recover the 
bulk STO nature. Although, similarly in the Ist-layer, 
U2xy,2xy is the largest on-site screened Coulomb inter- 
action in the 2nd-layer, the difference between U2xy.2xy 
and the others of t2g is only about 1 percent. Because the 
electrons confined at the interface have low dimensional- 
ity, thc correlation becomes effectively stronger than the 
bulk STO, even if the transfer integrals at the interface 
are similar to those of the bulk STO. 



4. Summary 

In this paper, we have determined the parameters of 
the low-energy effective model of LaAlOs/SrTiOa by 
the MACE. As with many interfaces, LaAlOs/SrTiOa 
has a complex band structure where the bands of both 
the interface and bulk regions are highly entangled with 
each other. By the disentangling scheme using thc max- 
imally localized Wannier function, we disentangle the 
low-energy part having the strong characters of Ti t2g 
orbital from the global KS-bands, and thus enable to 
calculate the parameters in the effective Hamiltonian of 
LaAlOa/SrTiOa by the cRPA. The parameters in this 
study offer not only the superlattice model but also 
the semi-infinite interface model (vacuum)-(LaA103)7v- 
(SrTi03)oo, because the screening effect from the inter- 
superccU rapidly decreases with distance and thus the 
superlattice model corresponds to the semi-infinite in- 
terface model in the limit of thick LAO and STO. The 
parameters have anisotropics and a layer dependence. 
The on-site energies in the Ist-layer are 650 meV higher 
than those in the 2nd-layer, which causes localization of 
the carriers at the interface. In the Ist-layer, while the 
transfer integral of the t2g orbital is similar to that of 
the bulk SrTiOs, the screened Coulomb interaction U of 
the t2g orbital becomes 10 percent larger than that of 
the bulk SrTiOg. In the bulk region of LaAlOg/SrTiOs, 
the parameters of the t2g orbitals recover the values 
and the symmetry of the bulk SrTiOs. The obtained 
3-dimensional parameters constitute the low-energy ef- 
fective models of LaAlOs/SrTiOa, either in semi-infinite 
structures in one of the directions or in supercell struc- 
tures. The resultant low-energy effective model offers a 
firm and quantitative basis when one wishes to solve the 
effective model by using accurate low-energy solvers in 
the future. 

Recently, an ab-initio dimensional downfolding 
scheme, which downfolds a 3-dimensional model to 
a lower-dimensional model in real space, has been 
formulated. '^^ By applying the dimensional downfolding 
scheme to the 3-dimensional model of LAO/STO, we can 
also obtain a 2-dimensional model of a single LAO/STO 
supercell. In the dimensional downfolding, one is able 
to expect that the weakness of inter-layer (inter-chain) 
couplings justify the RPA type perturbative treat- 
ment. In the present case, however, the interface layer 
is not particularly weakly coupled with the other layers. 
Nevertheless, the layers far away from the interface 
do not join in the low-energy excitations if the bulk is 
insulating as in the present case. Therefore, the gapped 
excitations in the bulk part can be safely downfolded 
into the metallic low-energy excitation near the interface 
in the same spirit of the cRPA. This enables to derive 
effective low-energy models solely for the interface part. 
Of course, we may need to keep the bands of not one 
but several layers near the interface in the low-energy 
models. This dimensional downfolding is a challenging 
future issue. 
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Table I. Transfer integrals for the t2g orbitals of the Ti sites in the bulk STO, tmn{Rx, Ry, Rz), where m and n denote symmetries of 
t2g orbitals. Units arc given in meV. 







(m, 


[0,0,0] 


[1,0,0] 


[1,1,0] 


(xy, xy) 


2430 


-299 


-109 


{^y, yz) 











{xy,zx) 











{y^,yz) 


2430 


-39 


7 


{yz,zx) 








9 


{zx, zx) 


2430 


-299 


7 



Table II. Transfer integrals for t2g orbitals of Ti sites in LA01.5ST03.5 calculated in LDA and GWA level, tmn{Rx, Ry, Rz), where m 
and n denote layers and symmetries of t2g orbitals. Units are given in meV. The parameters calculated from the LDA band are listed 
as "LDA" , and the parameters calculated from the LDA band with the self-energy correction of the La 4/ are listed as "GW" . 



LAO/STO 


LDA 


GW 


(m, 


[0,0,0] 


[1,0,0] 


[1,1,0] 


[0,0,0] 


[1,0,0] 


[1,1,0] 


{Ixy, Ixy) 


2606 


-289 


-101 


2607 


-288 


-104 


{Ixy, lyz) 





-7 


-1 





-4 


2 


{Ixy, Izx) 








-1 








2 


{Ixy, 2xy) 


-42 


6 


-9 


-42 


6 


-9 


{lxy,2yz) 





10 


6 





9 


6 


(Ixy, 2zx) 








6 








6 


{lyz, lyz) 


2617 


-57 


-2 


2634 


-59 


-5 


{lyz, Izx) 








-6 








-9 


{lyz, 2xy) 





8 


8 





8 


8 


{lyz, 2yz) 


-296 


3 


-6 


-299 


4 


-5 


{lyz, 2zx) 








8 








8 


{Izx, Izx) 


2617 


-302 


-2 


2634 


-298 


-5 


{Izx, 2xy) 








8 








8 


{Izx, 2yz) 








8 








8 


{Izx, 2zx) 


-296 


-98 


-6 


-299 


-98 


-5 


{2xy, 2xy) 


3268 


-298 


-109 


3264 


-298 


-109 


{2xy, 2yz) 





-3 








-3 





{2xy, 2zx) 




















{2yz, 2yz) 


3262 


-39 


6 


3262 


-39 


6 


{2yz,2zx) 








8 








7 


{2zx,2zx) 


3262 


-297 


6 


3262 


-297 


6 
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Appendix: GWA for disentangled band 

The quasiparticle energies and wave functions are ob- 
tained by solving 

{T + V,^t + Vu)^nk{r)+ J dr'E{r,r';Enk)Mr') 

(A-1) 

where T is the kinetic energy operator, and 14xt is the 
external potential, Vh is the Hartree potential. If we con- 
sider only the diagonal parts of the self-energy S, eq. 
(A-1) is reduced to 

Enk = ej;^^ - (?ife|yL°A|nfc) + (nfc|E(£;„fc)|nfc). (A-2) 

The self-energy operator must be estimated at the quasi- 
particle energy Enk- This is done by expanding the ma- 
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Tabic III. Effective Coulomb interaction (t/)/exchange (J) interactions between the two electrons for all the combinations of Ti-t2g 
orbitals in the bulk STO and LA01.5ST03.5, respectively (in eV). Especially, for the combinations between the orbitals in the same 
Ti sites, the matrix elements represent the effective on-site Coulomb interaction. The band structures of LA01.5ST03.5 are calculated 
in the LDA and GWA. 



STO 








U 














J 








xy 


yz 


zx 










xy 


yz 


zx 








xy 


3.76 


2.81 


2.81 








xy 




0.46 


0.46 








yz 


2.81 


3.76 


2.81 








yz 


0.46 




0.46 








zx 


2.81 


2.81 


3.76 








zx 


0.46 


0.46 










LAO/STO(LDA) 








u 














J 








Ixy 


lyz 


12:2; 


2xy 


2yz 


2zx 




Ixy 


lyz 


Izx 


2xy 


2yz 


22a; 


Ixy 


3.65 


2.63 


2.63 


0.49 


0.55 


0.55 


Ixy 




0.46 


0.46 


0.00 


0.00 


0.00 


lyz 


2.63 


3.48 


2.58 


0.55 


0.65 


0.63 


lyz 


0.46 




0.44 


0.00 


0.01 


0.00 


Izx 


2.63 


2.58 


3.48 


0.55 


0.63 


0.65 


Izx 


0.46 


0.44 




0.00 


0.00 


0.01 


2xy 


0.49 


0.55 


0.55 


3.62 


2.67 


2.67 


2xy 


0.00 


0.00 


0.00 




0.46 


0.46 


2yz 


0.55 


0.65 


0.63 


2.67 


3.59 


2.65 


2yz 


0.00 


0.01 


0.00 


0.46 




0.46 


2zx 


0.55 


0.63 


0.65 


2.67 


2.65 


3.59 


2zx 


0.00 


0.00 


0.01 


0.46 


0.46 




LAO/STO(GW) 








U 














J 








Ixy 


lyz 


Izx 


2xy 


2yz 


22a; 




Ixy 


lyz 


Izx 


2xy 


2yz 


22a; 


Ixy 


4.00 


2.98 


2.98 


0.72 


0.79 


0.79 


Ixy 




0.46 


0.46 


0.00 


0.00 


0.00 


lyz 


2.98 


3.83 


2.93 


0.78 


0.89 


0.87 


lyz 


0.46 




0.44 


0.00 


0.01 


0.00 


Izx 


2.98 


2.93 


3.83 


0.78 


0.87 


0.89 


Izx 


0.46 


0.44 




0.00 


0.00 


0.01 


2xy 


0.72 


0.78 


0.78 


3.80 


2.84 


2.84 


2xy 


0.00 


0.00 


0.00 




0.46 


0.46 


2yz 


0.79 


0.89 


0.87 


2.84 


3.76 


2.83 


2yz 


0.00 


0.01 


0.00 


0.46 




0.46 


2zx 


0.79 


0.87 


0.89 


2.84 


2.83 


3.76 


2zx 


0.00 


0.00 


0.01 


0.46 


0.46 
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Fig. A-1. Real and imaginary parts of the matrix elements of the 
self-energy S for a 4/ Wannier band in LA01.5ST03.5 at F- 
points. Solid and dashed lines represent the real and imaginary 
parts of the S, respectively. 



trix elements of the self-energy operator to the first order 
in the energy around e^k- Then the quasiparticlc energy 
is obtained explicitly; 

Enk = fnfc + ^rifc(erifc)(ej^fc^ — Ejifc + AE„fe(e„fe)), (A-3) 



where the sclf-encrgy is given by 



nh 



LDA 



and is the renormalization factor 



^-"^^^ ) ■ 



(A-4) 



(A-5) 



If e„fc coincides with e^^^, then eq. (A-3) is simplified; 

We show the real and imaginary parts of the self- 
energy for a 4/ Wannier band calculated by the 1-shot 
GW approximation in Fig. A-1. The real and imaginary 
parts of the self-energy strongly oscillate above the Fermi 
level Ep = eV, especially at 2-3 eV near the La 4/ 
bands in the LDA level. This lack of smoothness of the 
self-energy S is obtained in the "1-shot" correction. The 
self-energy S should hopefully be calculated in a self- 
consistent procedure of the Hedin's set of coupled equa- 
tions.'^^ In this calculation, we have calculated the self- 
energy correction by one shot without iteration to save 
the computational cost. Considering the unstable behav- 
ior of the self-energy near the La 4/ bands, we expand eq. 
(A-2) around Ep = eV and then quasiparticle energies 
are approximated in the first order as follows; 



(A-7) 



The approximation eq. (A-6) is equivalent to eq. (A-7), 
which is justified if the linearity of self-energy E in 
the low-energy region is eventually recovered in the 
self-consistent accurate estimate. The smoothness of E 
should be eventually obtained after the self-consistent 
calculation of the GWA. In this calculation, we average 
^nfc(O) and AS„fc(0) in terms of the band index and the 



fc-points; 

Enk = {Znk{0))nk{^n^^ 



(AS„fc(0))„i 



(A-8) 



Since wc are interested in the La 4f level measured 
from the Fermi level, the self-energy correction to the 
La 4f level is corrected by subtracting the correction to 
the Fermi level. The latter is approximately evaluated by 
the GW self-energy for the bottom of the t2g band us- 
ing eq.(A-7). Figure 6 shows the band structure thus ob- 
tained after the 1-shot GW corrections for the 4/ bands. 
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